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ABSTRACT 


A  method  is  presented  for  determining  the  location  of  an  uncooper¬ 
ative  radar  transmitter  by  timing  the  pulses  of  the  transmitter  after 
they  have  been  reflected  or  transduced  by  a  moving  target  of  known  loca¬ 
tion.  A  procedure  is  discussed  for  assessing  the  errors  involved. 

1.  THE  PROBLEM 

Ttie  following  problem  has  been  proposed  in  the  radiolocation  field; 
Given  a  radar  transmitter  in  a  foreign  location  tracking  a  missile, 
usually  a  missile  of  foreign  origin,  assume  that  at  a  receiving  point 
in  a  friendly  location,  it  is  possible  to  detect  the  radar  pulses  re¬ 
flected  from,  or  transponded  by,  the  missile  but  it  is  not  possible 
to  hear  the  transmitter  directly.  The  flight  path  of  the  missile  is 
known  at  the  receiver  location,  usually  through  friendly  optical  or 
radar  tracking.  The  problem  is  to  locate  the  radar  transmitter. 

Several  assumptions  must  be  made  in  order  to  solve  this  problem; 

The  pulse  repetition  period  T  of  the  radar  is  taken  as  constant,  although 
its  value  need  not  be  known.  If  the  pulse  is  transponded  at  the  missile, 
the  time  delay  of  the  transponder  is  also  taken  to  be  constant. 

The  location  of  the  transmitter  will  be  deduced  from  the  missile 
position  Information  together  with  data  giving  the  time  of  arrival  of 
the  reflected  or  transponded  pulses  at  the  receiver. 

2.  THEORETICAL  DISCUSSION  AND  SOLUTION 

Consider  the  times  of  arrival  T  of  the  received  pulses.  The 
receiver  time  at  which  the  nth  pulse  left  the  transmitter  (fig.  1)  is: 


total  transmit ter- to-mis si le-to-receiver4dl stance 

velocity  of  light 

time  delay  of  transponder 

(6=0  for  a  reflected  signal) 


where  H  = 
n 

c  = 
6  = 


Naturally,  the  missile  position  used  in  calculating  H  is  that  position 
the  missile  occupies  when  the  nth  pulse  is  operated  upon. 

Consider  another  pulse,  the  jth  pulse;  it  is  clear  that  it  left  the 
transmitter  a  period  of  time  (n-j)  T  before  the  nth  pulse,  provided  n  >  j. 
Therefore 
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Figure  I.  Tiae  disposition  of  transaltted  pulses 


=  m-i  . 


(n-j)  T 


Now  take 
It  follows  at 


two  other  pulses, 
once  that 


m  and  1,  m  >1,  such  that  n-j 


^'^n  ■  “  'fir 


T,)  =  H 


H.  -  H 
J  m 


Each  H,  say  H  ,  consists  of  two  parts:  a  distance  L  from  the 
receiver  to  the  missile,  which  is  known  to  the  receiver,  and  a  dis¬ 
tance  r  from  the  transmitter  to  the  missile,  which  is  not  known  at 
n 

the  receiver. 


H  = 
n 


+  L 


Using  this,  we  can  write 


(t  -  X.  -  X 
n  j  m 


■ 


+  L. 
J 


L  * 
m 


r ,  - 
J 


+  r. 


This  equation  will  be  called  the  trace  equation,  and  surfaces  or  curves 
in  space  that  satisfy  it  will  be  called  traces.  In  this  form,  the  left 
hand  side  of  the  trace  equation  consists  entirely  of  known  quantities. 
Note  that  we  may  choose,  if  we  wish,  j  =  m. 

Before  further  discussion,  a  coordinate  system  will  have  to  be 
introduced.  Although  one  might  be  tempted  to  choose  spherical  coor¬ 
dinates  because  this  problem  takes  place  on  the  earth's  surface, 
the  unknown  r's  are  distances,  and  the  expression  for  distance  is  gen¬ 
erally  simpler  in  rectangular  coordinates.  Therefore,  three-dimensional 
rectangular  coordinates  (X,  Y,  Z)  will  be  used. 

If  the  location  of  the  transmitter  is  known  approximately,  it  is 
advantageous  to  orient  the  coordinate  system  so  that  the  Z  coordinate 
of  the  transmitter  is  reduced  to  zero  or  very  nearly  so.  That  is,  the 
XY  plane  should,  cut  the  surface  of  the  earth  near  an  expected  location 
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of  the  transmitter.  This  procedure  has  the  advantage  of  reducing  the 
problem  to  a  two-dimensional  one,  which  is  more  tractable,  although 
this  simplification  is  not  essential. 


A  useful  suggestion,  due  to  Mr.  J.  E.  Tompkins  is  first  to  calcu¬ 
late  one  or  two  sets  of  data  on  a  three-dimensional  basis,  using  an 
arbitrary  orientation  of  coordinates.  These  results  may  then  be  used 
to  reorient  the  coordinate  system  so  as  to  reduce  the  problem  to  a  two- 
dimensional  one. 

Let  the  receiver  be  at  the  origin  (0,0,0);  the  unknown  transmitter 
location  will  be  (a,  3,  y)  •  ’rhe  numbers  (X  ,  Y  ,  Z  )  will  denote  the 
position  of  the  missile  when  it  operates  upon  t8e  n?h  pulse.  With 
these  conventions,  the  trace  equation  becomes 


-  7(x-X,)^  +  (y-V,)“  ♦  -7(x-X  +  (y-Y  )*  +  (z-z 

J  J  J  m  m  m 

+  +  (y-Y^)^  +  (z-Z^)^ 


The  trace  surfaces  (or,  in  degenerate  cases,  curves)  always  con¬ 
tain  the  point  (a,  3,  y)>  any  choice  of  n,  j,  m,  1,  subject  to 
the  condition  n-j  s=  m-i.  Since  the  surfaces  are  generally  not  identi¬ 
cal  parallel  or  nested  for  different  choices  of  the  subscripts,  then 
they  must  intersect.  In  the  case  of  surfaces,  at  least  three  must  be 
used  to  define  a  point  In  case  a  few  surfaces  (or  curves)  have  mul- 
tiple^lntersections ,  then  more  surfaces  must  be  drawn  to  eliminate  the 
spurious  intersections  from  the  solution.  Ehccept  probably  in  the  case 
where  the  missile  path  is  a  straight  line,  the  point  (a,  3,  Y)  will 
normally  be  the  only  point  of  intersection  of  all  the  surfaces  or 
curves. 

In  brief,  then,  the  transmitter  is  located  by  finding  the  simul¬ 
taneous  solution  of  a  set  of  trace  equations.  If  the  physical  situa¬ 
tion  is  such  that  the  problem  can  be  studied  well  enough  in  two  di¬ 
mensions,  all  that  is  required  is  to  graph  the  real  solutions  of  each 
trace  equation,  and  select  the  common  intersection  point.  For- a  three- 
dimensional  problem,  solid  or  matrix  graphs  are  still  a  possibility, 
but  a  better  approach  might  be  to  program  the  set  of  trace  equations 
on  a  digital  computer. 

3 .  ERROR  TREATMENT 


It  is  desirable  to  estimate  the  error  that  will  occur  when  applying 
this  method.  Such  error  will  arise  from  imprecisions  in  the  measurements 
of  times,  coordinates,  or  distances  involved.  Let  O  be  any  one  of  the 


numerical  parameters  involved  in  a  given  trace  equation,  and  suppose 
that  an  estimate  of  the  error  of  measurement  of  <7  shows  that  a  lies 
within  the  bounds  K  <  a <  \.  The  symbol  0  will  denote  the  mean  value 
(in  any  sense)  of  these  measurements  of  the  parameter. 

In  this  procedure,  the  graphical  presentation  is  essential. 

First )  gfraph,  for  each  trace  equation,  the  surface  or  curve  using  the 
value  ff";  this  will  be  called  the  median  surface  (curve).  Now  as  a 
varies  from  0,  the  trace  surface  deforms  continuously  in  at  least  the 
6-space  in  which  x,  y,  and  z  are  complex  quantities.  The  new  surfaces 
created  by  varying  0  from  "c  will  be  called  the  varied  surfaces  (curves). 
This  smooth  variation  la  due  to  the  fact  that  the  trace  equations  are 
continuous  functions  of  all  variables  and  parameters  Involved.  For 
values  of  (x,y,z)  in  the  immediate  neighborhood  of  (a,  p,  y) ,  the  solu¬ 
tion  of  each  trace  equation  will  be  wholly  real. 

Consider  any  point  P  on  the  median  surface  (curve);  the  varied 
surfaces  (curves)  will  generally  lie  on  both  sides  of  the  median  curve, 
depending  on  which  way  a  varies  from  5.  Denote  the  distance  at  P  to  a 
varied  curve  by 


[P  -  p] 

and  arbitrarily  denote  the  distance  on  one  side  of  the  median  curve  +, 
on  the  other  side  -.  The  point  pwill  be  the  point  on  the  varied  sur¬ 
face  closest  to  P. 

For  each  P,  there  are  points  qCp,  r^p  on  the  varied  surfaces, 
such  that 

[P-q]  e  sup  [P-p] 

[P-w]  =  inf  [P-p] 


where  w  and  q  are  found  by  varying  o  within  the  prescribed  limits.  A 
surface  is  then  defined  by  the  q's  and  another  by  the  w's;  these  sur¬ 
faces  (or  curves)  are  called  extreme  surfaces.  The  volume  (or  area) 
they  enclose  will  be  called  the  error  volume  (area)  E,  ;  the  subscript 
denotes  that  it  was  generated  from  the  kth  trace  equation. 

A  little  thought  will  now  show  that  the  target  will  lie  in  the 
volume  (area)  A  given  by 

A  =  E,  nE.,D  .  .  .  (^E 
12  w 

if  there  were  w  trace  equations  employed  altogether. 
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If,  Instead  of  only  one  parameter,  several  parameters  have  error 
tolerances,  the  principle  remains  the  same.  The  parameters  are  varied, 
singly  or  In  any  combination,  within  the  limits  of  their  respective 
errors,  until,  a  q  and  w  are  defined  for  each  P,  and  the  extreme  surfaces 
constructed  as  befora. 

It  will  be  recognized  that  the  extreme  surfaces  do  not  necessarily 
coincide  with  any  of  the  varied  surfaces.  The  extreme  surfaces  are 
rather  to  be  considered  as  envelopes;  as  such,  they  must  satisfy  the  dif¬ 
ferential  equation  of  a  family  of  curves,  at  least  as  a  singular  solution 
The  envelope  need  not  be  a  member  of  the  family  Itself. 

In  concluding  the  general  discussion  of  the  problem,  attention 
should  be  called  to  the  assumption  that  the  pulse  repetition  frequency 
(prf)  is  constant.  Although  this  point  has  not  been  rigidly  investigated 
it  seems  reasonable  that  this  requirement  can  be  relaxed  somewhat.  For 
example,  if  the  prf  Jitters  around  an  average  value,  this  will  seemingly 
not  lead  to  any  great  error  over  a  large  number  of  observations .  It  is 
also  very  likely  that,  if  data  be  examined  over  any  long  time  interval 
during  which  the  prf  has  a  time  average,  useful  results  will  be  obtained, 
even  if  the  momentary  excursions  of  the  prf  are  large.  Even  if  some  er¬ 
ror  is  caused  by  nonconstancy  of  the  prf,  in  some  applications  even  im¬ 
precise  results  are  better  than  no  information  at  all. 
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